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ABSTRACT 

Aims. We report on the results of an 8-year survey carried out at the La Silla Observatory with the HARPS spectrograph to detect and 
characterize planets in the super-Earth and Neptune mass regime. 

Methods. The size of our star sample and the precision achieved with HARPS have allowed the detection of a sufficiently large 
number of low-mass planets to study the statistical properties of their orbital elements, the correlation of the host-star metallicity with 
the planet masses, as well as the occurrence rate of planetary systems around solar-type stars. 

Results. A robust estimate of the frequency of systems shows that more than 50% of solar-type stars harbor at least one planet of 
any mass and with period up to 100 days. Different properties are observed for the population of planets less massive than about 
30 M e compared to the population of gaseous giant planets. The mass distribution of Super-Earths and Neptune-mass planets (SEN) 
is strongly increasing between 30 and 15 M e . The SEN occurence rate does not exhibit a preference for metal rich stars. Most of the 
SEN planets belong to multi-planetary systems. The orbital eccentricities of the SEN planets seems limited to 0.45. At the opposite, 
the occurence rate of gaseous giant planets is growing with the logarithm of the period, and is strongly increasing with the host-star 
metallicity. About 14% of solar-type stars have a planetary companion more massive than 50 M e on an orbit with a period shorter 
than 10 years. Orbital eccentricities of giant planets are observed up to 0.9 and beyond. 

Conclusions. The precision of HARPS-type spectrographs opens the possibility to detect planets in the habitable zone of solar-type 
stars. Identification of a significant number of super-Earths orbiting solar-type of the Sun vicinity is achieved by Doppler spectroscopy. 
41 newly discovered planets with HARPS are announced in the Appendix of this paper, among which 16 Super-Earths. 

Key words. Stars: late-type - Planetary systems - Techniques: radial velocities - Techniques: spectroscopy - Methods: statistical 
analysis 



1. Super-Earths and Neptune-mass planets around 
solar-type stars 

Thanks to the increasing sensitivity of instruments optimized for 
the measurement of stellar radial velocities, planets with masses 
smaller than 2 M e have been detected. Over the past few years 
an impressive progress has been made for the detection of close- 
in super-Earths (planets with masses between 1 and 10 M e ) and 
Neptune-mass planets. At present, the number of low-mass plan- 
ets discovered is sufficient to study the statistical properties of 
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this rich sub-population and to estimate their occurrence fre- 
quency around solar-type stars. 

The first hint to low-mass planets orbiting solar-type stars 
on tight orbits was provided by the discoveries of the short pe- 
riod companions hosted by 55 Cnc ( |McArthur et al.||2Q04| > and 
H Ara ( |Santos et al.|2004a) . Both discoveries were the result of 
intensive and special- interest Doppler monitoring: the study of 
the complex dynamics for the first system and an asteroseismol- 
ogy campaign for the second. Orbital solutions and masses for 
these two low-mass planets have been subsequently revised as 
more measurements were gathered and new approaches applied 
for the data analysis. For 55 Cnc e, |Dawson & Fa brycky (2010) 
analyzed the impact of observation aliases on the solution de- 
termination and proposed a revised value for the period passing 
from 2.7 down to P — 0.736 days corresponding to a new mini- 
mum mass ni2 sin i — 8.3 M ffi ) for the planet. In the case of /j Ara 
[Pepe et aL (2007 1 revised the mass of the planet at P — 9.6 days 
to m.2 sin i = 10.5 M ffl in the context of a new orbital solution 
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implying a total of four planets. These first detections induced 
important changes in our observing strategy, leading now to the 
discovery of a wealth of low-mass planets with short orbital pe- 
riods. Measurements were made long enough to diminish the 
stellar acoustic noise, as well as series of measurements over 
consecutive nights to lower granulation noise, were the key for 
demonstrating the existence of an extremely rich population of 
low-mass planets on tight orbits, these planets being frequently 
members of multi-planetary systems. In this context, a few land- 
mark contributions of HARPS to the field can be recalled: 

- In 2006, the discovery of three Neptune-mass planets hosted 
by HD 69830 ( |Lovis et aL|2006) , an interesting K star sur- 
rounded by a prominent disk of dust. 

- In 2009, the discovery of a planetary system with three super 
Earths orbiting HD 40307 ( [Mayor et al.|2009b) . 

- In 2011, the characterization of a system with 7 planets on 
tightly-packed orbits around HD 10180 ( |Lovis et al.|2011b i. 
This is the planetary system with the largest number of 
planets known so far, with the exception of our own Solar 
System. 

- After 4 years of the HARPS survey ( |Mayor et al.|2003| >, the 
number of low-mass planets detected was large enough to 
allow for a preliminary analysis of the statistical properties 
of this specific population ( |Udry & San tos 2007 [ [Mayor & 



illustrated for exampl e by the four-planet s ystem around GJ 581p] 
([Bonfils et al |[2005| |Udry et al.|[2007| |Mayor et al.]|2009a| 



Udry 2008 



The mass distribution of small-mass planets was 
found to be rising towards smaller masses down to 10M ffi , 
where the measurements start to be dominated by detection 
limits. A bi-modal mass distribution is observed for low- 
mass planets and giants. Concerning the planet-star separa- 
tion, the pile-up around 3 days observed for hot Jupiters was 



not observed for the low-mass population (Udry & Santos 
|2007] >. Already with the very first detections of low-mass 
planets, it appeared that these planets were hosted by stars 



without overabundance of metals ( |Udry et al.|2006 1, in heavy 
contrast with the strong correlation between host-star metal- 
licity and occurrence frequency observed for gaseous giants 
( [Santos et aLl[200T] [2004bl |Fischer & Valenti||2005]irThis 
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specific point was revisited and strengthened by 
(2008] who carried out a homogenous metallicity determi 
nation of the HARPS sample stars. 

A first estimate of the occurrence rate of low-mass planets 
(e.g. below 30 M e , with periods shorter than 50 days) was 
proposed by Lovis et aL] ( 2009} . The estimated frequency 
was ~ 30%. Note that, based on the 77 Earth project car- 
ried out at Keck, a significantly lower estimate for the oc- 
currence rate of low-mass planets was proposed by |Howard] 
et al . I (|20 1 0J . Their estimate was of ~ 18 % for planets be- 



tween 3 and 30 M® and periods shorter than 50 days. Their 
analysis included a correction for detection biases. 

Recently, the Kepler space mission contributed to signifi- 
cantly to our knowledge of low-mass planet properties by pub- 
lishing a list of 1235 transit detections of small-size planet can- 
didates (Boruc kT et al.|20lT| l. In addition to the confirmation of 
the existence of the rich population of low-mass planets on tight 
orbits, a fantastic result of the mission resides in the high fre- 
quency of multi-transiting systems in extremely co-planar con- 
figurations. As for the HARPS survey, the deduce that a large 
number of stars host several planets. A stunning example of such 
systems is given by the 6 transiting planets around Kepler- 11 
( |Lissauer et al.|20lT| >. 

The present study mainly discusses the properties of the pop- 
ulation of low-mass planets hosted by solar-type stars. We note, 
however, that this sub-population also exists around M stars, as 



Forveille & al. 2011|l.For a discussion of the HARPS companion 



subprogram devoted to M stars we refer to Bon fils & al.| ( |201 1) . 

The HARPS and CORALIE surveys used in the present 
study are described in Sect. 2, with a discussion of the stel- 
lar samples, the corresponding radial-velocity precision, and the 
correction of the detection biases. The planetary systems at the 
basis of our discussion will be presented in Sect. 3 and the de- 
rived distributions of planetary masses, periods, and eccentric- 
ities in Sect. 4. The emphasis of this study is on the results 
obtained for the population of super-Earths and Neptune-mass 
planets. However, combining our HARPS data with results of 
the CORALIE planet-search program allows us to derive sta- 
tistical distributions for the whole domain of planetary masses, 
from less than two Earth masses up to the most massive giant 
planets. This approach will in particular allow us to compare 
the derived statistical properties as a function of planetary mass. 
Clearly, with the radial-velocity technique, the detection of low- 
mass planets on longer periods is strongly affected by detection 
biases. These detection biases have been carefully modeled and 
corrected. This bias correction is of great importance when we 
estimate the occurrence of planets as a function of mass and pe- 
riod. This two-variable distribution is one of the most direct way 
to compare observations and Monte-Carlo based planet popu- 
lation synthesis as e.g. developed by Ida & Lin ( 2004| [2005 
2008a]b]l or Mordasini et al. ( 2009b|a| l. The physics of plan- 
etary formation dependents on the metallicity of the accretion 
disk. The large number of super-Earths and Neptune-mass plan- 
ets discovered in our survey provides a unique sample to derive 
the distribution of the host-star metallicity as a function of planet 
masses (Sect. 5). One of the most important goals of the field is 
certainly the estimation of the number of Earth-type planets ly- 
ing in the habitable zone of solar-type stars. Still more important 
would be the preparation of a list of stars with such planets at 
short distance to the Sun, preparing for follow-up studies aim- 
ing at detecting life signatures on those planets. High-precision 
Doppler surveys will contribute to the build up of such an "input 
catalogue" (see Sect. 6). 



2. The HARPS and CORALIE planet-search 
programmes 

2.1. Stellar samples 

2.1.1. The CORALIE survey 

A volume-limited radial -velocity survey is going on since 1998 
with the CORALIE spectrograph, located at the Nasmyth focus 
of the 1 .2-m EULER Swiss telescope at La the Silla Observatory. 
CORALIE is a twin of the ELODIE spectrograph (Ba ranne et al.| 
1996) running in the past at the Haute-Provence Observatory, 
and having allowed the discovery of 51 Pegb (Mayor & Queloz 
1995). In spite of the small telescope size, this instrument has 
a significant share of exoplanet discoveries. The volume-limited 
survey includes some 1650 stars of the southern sky. The limit 
distance for the stars in the sample is depending on the stel- 



Vogt et al. (20101 announced the detection of 2 additional planets 



in the G1581 system, with periods of 36.6 and 433 days, based on the 
analysis of the 119 published HARPS radial velocities (Mayo r et al.| 
2009a) and 122 Keck HIRES measur ements. With a final se t of 240 
high-precision HARPS measurements, Forveill e & al.| ( [201 1| > are now 
able to rule out the existence of the claimed 2 additional planets in the 
system. 
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lar spectral type: it has been set at 50parsecs for the F and G 
dwarfs, and linearly decreasing between KO and MO, in order 
to keep a reasonable exposure time to reach a given precision 
( |Udry et al.|2000| ). During the first ten years of observations, be- 
tween 1998 and 2008, the typical precision of CORALIE radial- 
velocity measurements was about 5ms _1 . In 2008, an upgrade 
of the spectrograph allowed for an improvement of its efficiency 
and precision. The goal of the CORALIE survey is to get, for an 
unbiased volume-limited sample, properties of substellar com- 
panions of solar-type stars, from sub-Saturn mass giant planets 
up to the binary star regime. Some interesting results have al- 
ready been obtained from this material as e.g. the correlation 
between host star-metallicity and occurrence frequency of giant 
planets (Santo s~et al.|200l" 2004b i, or a better definition of the 
location of the brown-dwarf desert, and hence of the observed 
maximum mass of giant planetary companions at intermediated 
periods (Sahl mann et al.|2010| l. 

2.1.2. The HARPS survey 

The HARPS instrument is a high-resolution spectrograph opti- 
mized for precise radial-velocity measurements, and specially 

2003). 



developed for exoplanet search programs (Mayor et al 



Temperature and pressure are controlled in a vacuum tank lo- 
cated in a thermal enclosure, the instrument is located in the 
gravity invariant Coude room and is fed with optical fibers is- 
sued from the Cassegrain adapter at the telescope. Thanks to its 
extreme thermal stability and the stable illumination of the op- 
tics provided by optical fibers (including a double-scrambling 
device), an instrumental precision better than 1 ms _1 is routinely 
achieved ( |Pepe & Lovis||2008| l. Both instruments (CORALIE 
and HARPS) are determining radial velocities by using the so- 
called cross-correlation technique (Baranne et al. 1996) between 



the extracted and calibrated high-resolution spectra and a binary 
(0,1) template. The combination of the instrument and the tech- 
nique is very efficient. For illustration, HARPS on the 3.6-m tele- 
scope reaches, on a 7.5 magnitude star, a precision of 1 ms _1 in 
only one minute. The dispersion relation for the orders of the 
stellar spectra on the CCD are determined by using a thorium- 
argon calibration lamp. A second parallel optical fiber is used 
for a simultaneous measurement of the calibration lamp during 
science exposures, monitoring thus potential drifts of the spec- 
trograph during the night, and allowing subsequently for a cor- 
rection of this drift if necessary. However, the nightly drift is 
typically much smaller than 1 ms 1 in HARPS. 

As a return for the construction of the instrument, 500 nights 
of guaranteed-time observations (GTO) have been granted by 
ESO to the HARPS scientific consortium in order to conduct a 
comprehensive search for exoplanets in the southern sky. As part 
of the different sub-programs carried out within the GTO since 
October 2003, a fraction of 50 % of the GTO time has been ded- 
icated to the exploration of the domain of very-low mass planets 
around solar-type stars (from late-F to late-K dwarfs), a domain 
of masses which did not yet present any detections in 2003. This 
program proved to be extremely successful and, after the end 
of the allocated GTO observing time (PI: M.Mayor) in April 
2009, it has been continued as an ESO Large Program (PI: S. 
Udry). The targets followed at high-precision with HARPS have 
been selected as a subsample of the CORALIE volume-limited 
sample, with additional constraints however: F and G stars with 
chromospheric activity index larger than logR' HK = -4.75 (or 
larger than -4.70 for K stars) or high-rotation star have been re- 
jected, as well as binaries. We were finally left with a sample of 
376 supposedly non-active stars. 



An additional exploratory program aiming at finding still 
smaller-mass planets around a small sub-sample of 10 stars in 
the GTO high-precision survey by applying a more demand- 
ing observing strategy to average the perturbing effects of stel- 
lar intrinsic phenomena (pulsations, granulation, activity) on the 
radial-velocity measurements (Dumusque et al. 201 lb\ is also 
being carried out on HARPS (PI: FPepe). This program has al- 
ready demonstrated the capability of HARPS to detect smaller- 
amplitude planets (Pepe et al. 201 1 1 as, e.g., a planet with radial- 
velocity amplitude as small as 0.56 ms" 1 and a planet in the hab- 
itable zone of a K dwarf (see also Sect. 6). 

The statistical results presented in the following are derived 
from the combined available data from these programs as per 
July 201 1. Altogether the analyzed data are based on more than 
450 observing nights with the HARPS spectrograph . 



2.1 .3. A volume-limited sample for the analysis 

As the data are delivered by two (or three, if we consider the 
Coralie Upgrade as a new instrument) different instruments with 
measurements of different precisions, as well as of different ob- 
servation time spans and strategies, we adopted a statistical treat- 
ment taking into account these various constraints and allowing 
us to use in an optimum way the information content of each of 
the surveys. The CORALIE survey, with its more modest radial- 
velocity precision but longer duration, mostly provides informa- 
tion about gaseous giant planets with orbital periods covering 
a long time interval from below 1 day to more than 12 years. 
The HARPS survey, on the other hand, provides a very deep in- 
sight into the domain of Super-Earths and Neptune-type planets 
up to periods of about 100 days. To avoid as much as possible 
differential biases between the 2 samples, the CORALIE sample 
used in the analysis has been re-defined following the same cri- 
teria for non-active star selection as for the HARPS sample. As 
stellar activity does not a priori correlates with distance to the 
Sun, the sample keeps its statistically well-defined characteris- 
tic. Finally, the global sample for the statistical analysis includes 
822 non-active stars allowing to estimate detection limits, ob- 
servational biases, and correct for them in a large domain in the 
W2sin/ - logP plane. The HR-diagram of the final sample is 
shown in Fig.[T] 



2.1 .4. Detection limits 

Since every star in the statistical sample has a different observ- 
ing history, we computed, as a function of period, the max- 
imum planet mass of a potential planet which still compati- 
ble with the RV measurements of this same star. For each star 
(with or without a detected planet candidate) the so computed 
detection limit determines if it usefully constrains the distribu- 
tion at a given mass and period. Several authors (Cummings & 
Deliyannis||2007t |Zechmeister et al.|2009| |Howard et al.||2010[ 



Bonfils & al.||2011| i have presented statistical analysis of their 
planet-search program and discussed methods for finding unbi- 
ased distributions consistent with a set of detections. They all 
assume circular orbits to estimate the upper limits on the ampli- 
tude. We also adopt this safer approach to avoid problems in case 
of inadequate phase coverage. Moreover, |Endl et aL] ([2002 ) have 
shown that the detectability is only slightly affected by this hy- 
pothesis for eccentricities below e = 0.5. The originality of our 
method is to combine for each star up to three data sets from the 
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Fig. 1. HR-diagram of Hipparcos 2008 catalogue (black) , 
CORALIE volume limited sample (red), HARPS sample of low 
activity stars (green) and stars with planetary systems (blue) 



HARPS and CORALIE spectrographs^] In order to constrain the 
radial velocity zero point, times series with at least four points 
on each instrument are combined. If the star is known to host a 
planetary system, the velocity offset is left as a free parameter 
of the fitting process. Otherwise a% 2 test is used to determine if 
a linear trend represents a significant improvement compared to 
the constant model. After combining the data sets and removing 
the identified planetary signals, the stellar activity correlation is 



removed using the relation with the log R' HK described in Lovis 



et al. (201 la i. We obtain this way a single set of residuals per 
star which has been cleaned from detected planetary signatures, 
identified activity cycles and instrumental offsets. A search for 
a residual excess of periodic variability is then carried for every 
star. 



We followed the procedure described in Zechmeister et al. 
(2009}: We used the generalized Lomb-Scargle (GLS) peri- 
odogram to explore FAP level above 10~ 3 over 5000 bootstrap 
randomization of the data set. Whenever a possible periodicity 
in the RV data is unveiled at an FAP of 1%, a careful search for 
correlation with stellar parameters (log R' HK , FWHM, line bisec- 
tor) is carried on. If the signal doesn't seem to be related to the 
stellar activity, we use its period as a guess for a Keplerian fit 
with Yorbit (Segransan et al.||20lT) . If no credible orbit can be 
found or if a correlation with stellar jitter is suspected, the signal 
is removed by fitting a sine wave. Otherwise, the planet orbit is 
subtracted and the signal flagged as possible planetary candidate 
to be analyzed further. The whole procedure is re-applied on the 
so-obtained new residuals until no more significant signal can be 
found. It must be note that planet candidates found during this 
procedure (a total of 6) are considered for the statistical study 
but are non listed in the table of detected planets (see Table[5]l. 



2 The efficiency of CORALIE was increased in 2007 through a hard- 
ware optimization that has affected the radial velocity zero point. For 
this reason, we consider the upgraded spectrograph as a new instrument. 



After this operation we are able to derive for every star the 
period-mass limit above which we can exclude or detect a planet. 
As all the significant signals where removed from the time series, 
the residuals can be considered as pure noise. On these residu- 
als we first apply the GLS on 5000 bootstrap randomizations 
of the data set to determine a power threshold corresponding to 
1% FAP. The presence of a planet at the trial period P is then 
mimicked by adding a sine wave to the data. The calendar of 
measurements remains unaffected. The planet is considered to 
be detected if the power in the GLS periodogram at the period 
P (or at its 1 day alias) is greater than the threshold. Using the 
corresponding sine wave amplitude and the stellar mass, we can 
deduce the mass corresponding to a 99% detection limit. Due 
to the imperfect sampling of the measurements, the detection 
limit depends in practice on the phase of the signal. Therefore 
we compute the detection limits over 20 equidistant phases be- 
tween and 1 and treat them as 20 different stars. This solution 
provides a simple way of 'averaging' over phases with variable 
detection limits. 



2.1 .5. Completeness of the survey and planetary rate 

The global efficiency of the survey, as represented in Fig 5]and|7] 
is deduced from the upper mass envelope of the independent lim- 
its of the sample. For example, the 20% line represents the limit 
below which only 20% of the stars have a 99% detection. Such a 
diagram displays the fraction C(m sin i, P) of stars with sufficient 
measurements to detect or exclude a planet at a given period and 
mass. In other terms it tells us what is the subsample of stars 
useful to estimate the planetary rate in any mass-period region 
and how many planetary companions were potentially missed 
due to insufficient constraints provided by the observations. In 
order to determine and compare the planetary occurrence rate in 
different regions of the diagram, we assign an effective number 
Njj — 1 / Cfmsini, P) for each detected planet i around the star j. 
The planetary rate is then simply given by f p i = -X- £ ; N/j where 
N„ is the total number of stars in the sample. The sum is calcu- 
lated over every detected planet i in a given mass-period domain. 
Computing the occurrence of planetary systems is similar, using 
fsyst = -fr TjjN max j, except for the sum which is calculated on 
every host star j accounting only for the companion with the 
largest effective number N/j in the trial domain. 



3. Planetary systems in the combined sample 

155 planets belonging to 102 planetary systems have been iden- 
tified among the stars in our combined sample. More than two 
third of these planets have been discovered by the CORALIE 
and HARPS surveys. HARPS or CORALIE orbits also exist for 
most of the planets in our sample which were discovered in the 
frame of other programs. The list of the 102 planetary systems 
with orbital parameters and references to the discovery papers is 
provided in Table[3] The characteristics of the newly discovered 
planets and planetary systems included in the statistical analy- 
sis are furthermore presented in the appendix of this paper. The 
detailed analysis of these systems will be presented however in 
the papers referred to in Table[3] For the statistical analysis, the 
6 additional 'planetary candidates' issued from the computation 
of the detection limits have been included. 

A key parameter for the statistical discussion of low-mass 
planets is the global efficiency of the Doppler measurements in a 
given survey. This efficiency includes the precision of the spec- 
trograph, sample star characteristics, and the observing strategy 
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(including measurement frequency, exposure time long enough 
to diminish the different intrinsic stellar noises, and photon 
noise). In Fig. [2] we plot the representative points of the radial- 
velocity semi-amplitudes K as a function of the logarithm of the 
period log P for all super-Earths orbiting solar-type stars. Planets 
detected in the frame of other surveys are identified by red sym- 
bols and those of our surveys by blue dots. This plot is a straight- 
forward illustration of the contribution of HARPS in this field. 
The number of the newly detected super-Earths is impressive, as 
well as the very small amplitudes presently detected by the in- 
strument. This diagram shows how important the measurement 
precision turns out to be for the correct estimation of the occur- 
rence of low-mass planets, which may be strongly biased if the 
low-mass planet detections delivered by HARPS were lacking. 

Figure[3] presents the radial-velocity raw dispersion mea- 
sured on the HARPS targets (only the histogram of rms less than 
10ms _1 is displayed). Superimposed, we have also illustrated 
the histogram of the rms for stars with detected planetary sys- 
tems. We should note that some planetary systems are detected 
around stars having quite modest raw radial-velocity scatter, in 
some cases smaller than 2ms~ 1 . Once again, the importance of 
the measurement precision turns out to be of fundamental impor- 
tance when exploring the domain of super-Earth masses. After 
fitting planetary orbit(s) to our set of measurements we derive 
the residuals around the orbital solutions (O - C). These residu- 
als include all potential sources of noise, as instrumental errors, 
photon noise, stellar intrinsic noise (acoustic oscillation, granu- 
lation, activity, magnetic cycle), as well as still undetected small- 
mass planets. Magnetic cycles, analogues to the Sun's 11-year 
cycle, can furthermore induce variations of the rad ial velocities 
over several years and amplitudes up to 10ms _1 (Lovis et al 



|201 la] >. This effect is observed also for some stars with rather 
low activity level. Conveniently however, several spectral signa- 
tures (log R' HK , width of the cross-correlation function (CCF) or 
spectral line and CCF bisectors) can and must be used to identify 
magnetic cycles and correct for their influences on the velocities 
( |Lovis et al.|201 la| lDumusq ue et al.|201 la| ). These spectral sig- 
natures can also be used to test velocity variations on shorter 
time scales (days or weeks) to detect the possible influence of 
spots or other anisotropic features on the stellar surface. These 
diagnostics are of particular importance when searching for low- 
amplitude velocity variations, typical of low-mass planets. All 
these tests have been applied to the HARPS and CORALIE 
planet-host stars and are described in the corresponding discov- 
ery papers. The histogram of the O - C (Fig. [4] for the HARPS 
planetary systems is certainly the most significant estimation of 
the global precision of the program. The mode of the O - C dis- 
tribution lies at 1.4ms" 1 but many planetary systems show O — C 
smaller than 1 ms _1 even with a large number of measurements 
spanning several years of observations. 



4. Statistical analysis of the sample 

4.1. The mass-period distribution 

The (ni2 sin z, log P) diagram for the planetary systems and can- 
didates detected in our survey is displayed in Fig. [5] Fig. ^illus- 
trates the same information, including lines of detection prob- 
ability at various levels of completeness for the global sample. 
Several remarks can be made: 

- Objects with masses of about 1 ms _1 at very short period, and 
about 5 ms _1 at one year define a lower envelope of the dis- 
played planets. The planets at this lower envelope are located 
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Fig. 2. Radial-velocity semi-amplitude K as a function of orbital 
period for super-Earths (M < 10 M e ) hosted by solar-type stars. 
HARPS detection are plotted as blue dots and objects from the 
literature in red symbols (circles for the southern sky and square 
for the northern sky). 



in a region of very low-detection probability. The discovery 
of 9 planets close to the 2 % detection probability line was 
only possible thanks to an extremely large number of radial- 
velocity measurements of these stars. The detection of these 
planets is already a strong hint to the high occurrence rate of 
low-mass super-Earths. 

The density of massive super-Earths seems to exhibit a 
"mass-period" relation. The extremum of the density appears 
to drift from about 6M ffl at 10 days to more than 10M ffl at 
100 days. We do not believe that observation biases can ex- 
plain the upper envelope of this feature relation, since heavy 
planets at sorter periods would be detected more easily (as 
indicated by the detection limits). Given the fact that the low- 
mass planets are all issued from the HARPS sample, it is 
interesting to consider the detection limits achieved for the 
376 HARPS-stars only measured. The (m,2 sin i, log P) distri- 
bution for the HARPS subsample only is displayed on Fig. [7] 
We note that the average detection probability of a Neptune- 
mass planet with a period of 100 days is still 60%. The 
strongly decreasing density of super-Earths and Neptune- 
mass planets (almost coincident with the red line correspond- 
ing to a detection of 60 percent) seems to be a highly signif- 
icant feature. 

We do not have Neptune-mass planets (and objects up to 
50M ffl ) with periods larger than about 1000 days, despite the 
still non-neglectable detection probability in that domain of 
mass-period. For the moment this feature is not clearly un- 
derstood. It could be the result of a detection bias related 
to the lack of measurements done for these low-amplitude 
and long-period objects, or simply due to the paucity of such 
planets beyond 1 or 2 AU. For the moment we do not con- 
sider this feature as being significant. 
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Fig. 3. Histogram of the radial- velocity dispersion for the stars of 
the HARPS sample (black line). Only the part of the histogram 
with rms < 10 ms" 1 is displayed. The red histogram represents 
the raw radial-velocity dispersion for stars with detected exo- 
planets. 



- An "upper mass-log P" relation is observed in the diagram, 
as already pointed out in (Udry et al. 2003| l. The diagram 
suggests that the upper mass for jovian planets on short pe- 
riods (10 days) is close to one Jupiter and rises to 15 Jupiter 
masses at orbital periods of 10 years. 

- Not only the upper mass of planets is period dependent but 
the frequency of gaseous giant planets is strongly increasing 
with the logarithm of the orbital period, a feature already 
noticed in several giant planet searches ( [Marcy et al.|2005| 
Udry & Santos 2007). Figure[8] illustrates the histogram of 
the masses of planets with m.2 sin i > 50 M e as a a function of 
orbital period. The observed distribution (continuous black 
line) as well as the bias-corrected distribution (dashed red 
line) illustrate the rising occurrence of gaseous planets with 
orbital periods. Setting the low-mass limit of the considered 
subsample to 100M ffl does not change the general shape of 
the distribution. 



4.2. Occurrence rate of planets as a function of planetary 
masses 

The estimation of the frequency of stars with a planet in a spe- 
cific domain of the m.2 sin i - log P parameter space requires the 
correction of observations from the detection incompleteness. 
Actually as explained in Sect. 2.1.4 it is more taking into ac- 
count the incompleteness determination in our occurrence esti- 
mate. In the case we consider a large zone in the 1112 sin i - log P 
plane, the presence of multi-planetary systems in the zone makes 
this operation ambiguous. For those systems, we have adopted 
the correction associated with the planet the most difficult to 
detect, as its higher statistical weight dominates the area. We 
present in Table[T]the occurrence rate for different domains of 
masses and periods. We do not try to do any extrapolation of our 



RVO-C [m/s] 

Fig. 4. Histogram of the O - C residuals after fitting the (multi- 
) Keplerian model to the velocity measurements of a star with 
planets. The O - C dispersions include all possible sources of 
noises (photon noise, instrumental errors and stellar jitter), as 
well as the signal from still undetected planets. Note that a sig- 
nificant fraction of the O - C lies below 1 ms" 1 . 



results out of the range of the observed parameters, as e.g. to- 
wards smaller masses or longer periods. For example, as we do 
not have any knowledge of the distribution of orbital periods of 
super-Earths above ~ 100 days, we consider an extrapolation out 
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Fig. 5. Plot of the 155 planets (dots) and 6 candidates (triangles) 
of the considered HARPS +CORALIE sample in the m2 sinj - 
log P plane . 
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Fig. 6. Same as Fig. [5] with detection probability curves super- 
imposed. These detection probabilities are valid for the whole 
sample of 822 stars. After correcting for the detection bias, the 
fraction of stars with at least one planet more massive than 
50M ffi and with a period smaller than 10 years is estimated to be 
14 ± 2 %. The red points represent the planets which have been 
used to compute the corrected occurrence rate in the box indi- 
cated by the dashed line. The planets lying outside the box or 
being part of a system already taken into account are excluded; 
they are shown in black. 



to the habitable zone around the stars as being of no significance 
the time being. This is also true for the distribution towards very 
low masses, at any given period. An exponential extrapolation of 
the mass distribution towards very low mass planets is not jus- 
tified if, for instance, a third, distinct population of Earth-mass 
planets is supposed, as suggested by Monte-Carlo simulations of 
planet population synthesis based on the core-accretion scenario 
(e.g. |Mordasini et al.|2009b|a) . 

Based on the 77-Earth survey carried out at Keck observa- 
tory, |Howard et al.| ( |2010| ) have derived an estimate of the oc- 
currence rate of low-mass planets on tight orbits (< 50 days) as 
well as their mass distribution. The comparison of our results 
in the same range of parameter space with the efa-Earth esti- 
mate is given in Table|2] The results are also shown graphically 
in Fig. 9] where the occurrence rate is indicated for the boxes 
defined by |Howard et al.| ( |2010| l. As a result of our much larger 
stellar sample, but also due to the better sensitivity of HARPS 
with regard to low-mass planets, the number of detected plan- 
ets is significantly larger than the one issued from the Keck sur- 
vey. This is specifically important for super-Earths and Neptune- 
mass planets (see Fig. [2] Our estimate for the occurrence rate 
of planetary systems is stunningly large. Table[T] lists the esti- 
mation of the frequency of planetary systems for different part 
of the m.2 sin i - log P plane. We should emphasize for example 
that 75 % of solar type stars have a planet with a period smaller 
than 10 years. This result is limited to the domain of detectable 
planets, and does not involve extrapolations out of that domain. 
Another interesting result is obtained for the low-mass planets 
on tight orbits: about half of solar-type stars is an host of this 
type of planets. 



Period [days] 

Fig. 7. Same as Fig^but only for the HARPS subsample. The 
occurrence rate of planetary systems in the limited region be- 
tween 3 and 100M e , and with P < 1 year, is 47 ± 7 %. Again, 
only the red dots and the blue triangles (candidates) have been 
considered for the computation of the occurrence rate. 
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Fig. 8. Histogram of the planet frequency for planets with masses 
(m.2 sin i > 50 M®). The occurrence rate for gaseous giant planets 
is strongly increasing with the logarithm of the period log P. 



Despite the significant size of our sample of 822 stars, the 
number of hot Jupiters is quite small (5 planets with P < 
1 1 days). The estimated frequency, 0.9 + 0.4% is compatible 
with previous estimates. (For a more complete discussion of the 
CORALIE-alone sample we refer to Marmier et al. (in prep). We 
know that the mass distribution of planets hosted by G dwarfs 
is different from its equivalent for planets hosted by M dwarfs 
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Table 1. Occurrence frequency of stars with at least one planet in the defined region. The results for various regions of the m.2 sin i ■ 
log P plane are given. 



Mass limits 


Period limit 


Planetary rate based on 


Planetary rate 


Comments 






published planets 


including candidates 




>50M e 


< 10 years 


13.9 ± 1.7 % 


13.9 ± 1.7% 


Gaseous giant planets 


> 100 M e 


< 10 years 


9.7 ± 1.3 % 


9.7 ± 1.3% 


Gaseous giant planets 


>50M e 


< 1 1 days 


0.89 ± 0.36 % 


0.89 ± 0.36 % 


Hot gaseous giant planets 


Any masses 


< 10 years 


65.2 ± 6.6 % 


75.1 ± 7.4% 


All "detectable" planets with P < 10 years 


Any masses 


< 100 days 


50.6 ± 7.4 % 


57.1 ± 8.0% 


At least 1 planet with P < 100 days 


Any masses 


< 100 days 


68.0 ± 11.7% 


68.9 ± 11.6% 


F and G stars only 


Any masses 


< 100 days 


41.1 ± 11.4% 


52.7 ± 13.2% 


K stars only 


<30 M e 


< 100 days 


47.9 ± 8.5 % 


54.1 ± 9.1 % 


Super-Earths and Neptune-mass planets on tight orbits 


<30M e 


< 50 days 


38.8 ±7.1% 


45.0 ±7.8% 


As defined in Lovis et al. (2009 1 



Table 2. Comparison of detected planets detected and occurrence rate of the rj Earth survey (Howard et al. 2010 1 and HARPS- 
CORALIE survey. The comparison is restricted to planets with orbital periods smaller than 50 days. A^i stands for the numbers of 
detected planets, N2 stands for the number of candidates, and N3 represents the estimated occurrence rate of planets in the given 
mass range. 



Mass range 


HARPS & CORALIE survey 


77 Earth survey 


Me 


Nb of planets 

Ni N 2 N 3 


Planetary rate 
[*] 


Nb of planets 

iVi N 2 N 3 


Planetary rate 
[%] 


3-10 
10-30 
30-100 
100-300 
300-1000 


19 2 48.5 
25 1 20.6 
5 1 4.6 
4 0.8 
2 


16.6 ±4.4 
11.1 ± 2.4 
1.17 ±0.52 
0.58 ± 0.29 
0.24 ±0.17 


5 3 10.2 
4 1 4.6 

2 
2 
2 


11.8 ±4.3 

6.5 ±3.0 

1.6 ± 1.2 
1.6 ± 1.2 
1.6 ± 1.2 




4.3. The mass distribution 



Period 



[days] 

Fig. 9. Same as Fig.|6] The dashed lines represent the boxes in 
which the occurrence rate is computed as defined by |Howard| 
|et al.| l |2"0l0) . An additional box is shown for masses between 1 
and 3 Mm. 



(Bonfils & al. 201 l| l. Despite the rather limited range of stel- 
lar masses in our sample, we have tried a comparison of the 
(m.2 sin i - log P) distribution for dwarf stars of spectral type F 
and G versus the distribution for K dwarfs. The observed differ- 
ence of planetary rate (for periods smaller than 100 days) for the 
two spectral types is however not significant (Table[T]i. 



On Fig 10 we have plotted the histogram of masses of the plan- 
ets detected in our sample. We observe a drastic decline of 
the observed mass distribution from about 15 to 30 M B . If we 
limit the range of orbital periods and only consider planets with 
P < 100 days (Fig.[TT|), a region where the detection bias are not 
too important for low-mass planets, we immediately observe the 
preponderant importance of the sub-population of super-Earths 
and Neptune-mass planets in that domain of periods. After cor- 
rection of detection biases (Fig. 12 1, we see even more clearly the 
importance of the population of low-mass planets on tight orbits, 
with a sharp decrease of the distribution between a few Earth 
masses and ~ 40 M e . We note that the planet population synthe- 
sis models by Mordasin i et al.| ( |2009b| l predicted such a mini- 
mum in the mass-distribution at precisely this mass range. They 
also pointed out that a radial-velocity measurement precision of 
about 1 ms"' was required in order to detect this minimum. In the 
framework of the core accretion model, this can be understood 
by the fact that this mass range corresponds to the runaway gas 
accretion phase during which planets acquire mass on very short 
timescales. Therefore, unless timing is such that the gaseous disk 
vanishes at this moment, forming planet transits quickly through 
this mass range and the probability to detect these types of plan- 
ets is reduced correspondingly. In Fig. [12] the importance of the 
correction of the detection biases below 20M ffl is only the re- 
flection of the present observing situation for which only a lim- 
ited fraction of the sample has benefited from the large enough 
number of HARPS measurements, required to detect small-mass 
objects. Part of this correction is also related to the growing im- 
portance of the sin / effect with decreasing masses. 
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Fig. 10. Observed mass histogram for the planets in the com- 
bined sample. Before any bias correction, we can already notice 
the importance of the sub-population of low-mass planets. We 
also remark a gap in the histogram between planets with masses 
above and below ~30M ffl . 



4.4. The period distribution of Super-Earth and 
Neptune-mass planets 

The observed distribution of orbital periods for planets less mas- 
sive than 30M ffl is illustrated in Fig 13 In Fig 14 the same dis- 
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Fig. 11. Same as Fig. 



10 



) but for planets with periods smaller than 
100 days. We see the dominance of low-mass planet with short 
orbital periods. 



M2sini [Earth Mass] 

Fig. 12. Histograms of planetary masses, comparing the ob- 
served histogram (black line) and the equivalent histogram after 
correction for the detection bias (red line). 



tribution is reproduced with a black histogram, to be compared 
with the histogram after correction for detection incompleteness 
(red histogram). In agreement with Kepler's preliminary find- 
ings ( |Borucki et al.|20lT) , the sub-population of low-mass planet 
appears mostly confined to tight orbits. The majority of these 
low-mass planets have periods shorter than 100 days. Low-mass 
planets on longer periods are of course more affected by detec- 
tion limits, this is however, at least partly, taken into account in 
our bias estimate and correction. We conclude that this feature 
must be real. 



4.5. Orbital eccentricities of Super-Earth and Neptune-type 
planets 

Figure[l5] displays the orbital eccentricities as a function of the 
planetary mass. We can remark the very large scatter of orbital 
eccentricities measured for gaseous giant planets, some of them 
having eccentricities as large as 0.93. Such very large eccentric- 
ities are not observed for planets with masses smaller than about 
30 M ffl for which the most extrem values are limited around 0.45. 
For low-mass planets the estimation of small orbital eccentricites 
of the best keplerian fit is biased. For the moment, the eccentric- 
ities below 0.2 (and small masses) have to be considered with 
caution . 



4.6. Fraction of multiplanetary systems with low mass 
planets 

For systems with planets less massive than 30 M ffl , the fraction 
of multi-planetary systems is extremely high. For the 24 con- 
cerned systems this fraction exceeds 70 %. It is tempting to have 
a rate of multi-planetary systems hosting at least one gaseous 
giant planets. Unfortunately, the optimum observing strategy 
needed to detect low-mass planets has not been applied to every 
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Fig. 13. Observed period distribution for low-mass planets Fig . 15 . Mass-eccentricity diagram for the planets in the com- 
(m 2 sin i < 30 M e ) bined samp i e . 

characteristics can be listed: HD 10180, HD 1 1964, HD 134060, 
HD 160691, HD 181433, HD 204313, and HD 215456. 

5. Host star metallicities as a function of planetary 
masses 

The occurrence rate of giant gaseous planets strongly correlates 
with the host star metallicity. Large unbiased studies have pro- 
vided a robust and well-defined relationship between the fre- 
quency of gaseous giant planets and the metallicity of their host 
star (e.g. |Santos et aLl|2TjoT| |2004b[ |Fischer & Valenti|[2005"i. 
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The lack of correlation between the overabundance of heavy e 
ements and the mass of the convective zone of the star (Santos 
et al.|2004b] l is seen as a strong argument in favour of an origin 
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Fig. 14. Observed period distributions for low-mass planets 
(m.2 sin i < 30 M e ) before (black histogram) and after (red his- 
togram) correction for the detection bias. Most of the low-mass 
planets are confined on short period orbits. The mode of the dis- 
tribution appears just over 40 days. 



stars with giant planets in past. Presently, we observe a multi- 
planetary rate of only 26 % for these planets. This point will 
have to be revisited with additional precise velocities. For the 
formation of planetary systems, the existence of systems with 
one gaseous giant planet with large period and a low-mass planet 
on a tight orbit could be of interest. A few systems having these 



of the planet-metallicity correlation linked with the primordial 
abundance of the molecular cloud. 

For lower-mass planets, already after the very few first de- 
tections ( |Udry et al. 2006) suggested the absence of correlation 
between the host star metallicity and the presence of low-mass 
planets. This first claim was later confirmed by (Sousa et al. 
2008) on a larger but still limited sample. The present analy- 
sis of a statistically well-defined, much larger sample of low- 
mass planets, offers the opportunity to have a much more ro- 
bust insight in the relation between planet occurrence frequency, 
host star metallicity, and planet mass. For the analysis, multi- 
planetary systems will be only characterized by their most mas- 
sive planet. In Fig. [16] we can compare the histograms of host- 
star metallicities for planets with masses smaller than 30 M ffi and 
its equivalent for gaseous giant planets. On the same plot, the 
metallicity distribution for the stars in the global combined sam- 
ple is illustrated for an enlightening comparison. It is however 
not clear whether we observe a discontinuity in the-host star 
metallicity distributions with an increasing planetary mass. On 
Fig 17 the ([Fe/H], mi sin;) diagram can be used to set a limit 
between the two regimes of host star metallicity (if such a limit 
proves to be meaningful!). On the diagram, such a limit can be 
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Fig. 16. Histograms of host star metallicities ([Fe/H]) for giant 
gaseous planets (black), for planets less massive than 30 M© 
(red), and for the global combined sample stars (blue). The lat- 
ter histogram has been multiplied by 0. 1 for visual comparison 
reason. 



set at about 30-40 M e . At the exception of a single star, all the 
stars hosting planets less massive than 40 M ffl have metallicities 
below [Fe/H] = 0.20. This is well in contrast with the situation 
for stars hosting more massive planets. Interestingly this mass 
corresponds about to local the minimum in the mass distribu- 
tion between Neptunes and gaseous giants (see e.g. 
Fig.[TT}. 

The correlation between the occurrence of giant planets and 
the metallicity of the host star (i.e. the metallicity of the ma- 
terial in the proto-planetary disc) is a natural outcome of the 
core accretion planet formation theory. In this paradigm, mas- 
sive gaseous planets form by runaway gas accretion onto cores 
exceeding a critical mass, typically of the order of 10 - 20M ffl . 
The gas accretion from the disc goes on until the disc vanishes, 
typically after a few million years. Hence, the sooner in the evo- 
lution of the disk a critical-mass core can form, the larger the 
amount of gas that will still be available for accretion. A high 
metallicity (interpreted as a large dust-to-gas ratio in the mod- 
els) and/or massive discs favors the early growth of such critical 
cores. Conversely, lower-mass planets that do not accrete signifi- 
cant amount of gas, can grow their cores over a longer timescale 
and therefore do not depend as critically upon the metallicity. 
These effects have been born out in the population synthesis 
models by |Mordasini & al.| ( [201 1) . 



6. Planets in the habitable zone of solar-type stars 

All the very specific properties of the population of low-mass 
planets (super-Earths and Neptune-mass planets) are of special 
interest for constraining the formation of planetary systems. In 
addition to this, the surveys targeting low-mass planets, in the 
same way as the parallel efforts aiming to increase the precision 
of spectrographs optimized for Doppler measurements are pur- 



M2sini [Earth Mass] 

Fig. 17. Estimation of the planetary-mass limit between the two 
regimes for the metallicity dependance of host stars. A vertical 
line at 30 M e distinguishes the two populations. We should note 
that such a limit also corresponds to the gap in the mass distribu- 
tion (see FigffO|and[T2]i. On the right side of the vertical line we 
do not observe significant changes of the metallicity distribution 
above 30 Mffi. We remark that stars with metallicity exceeding 
0.15 are for their huge majority associated with planets more 
massive than 30M ffl . 



suing a still more difficult challenge, the detection of Earth twins, 
rocky planets orbiting stars in the so-called habitable zone, and 
possibly around stars as close as possible to the Sun. This last 
condition is of special importance for future experiments aiming 
at the spectroscopic follow-up of the planet, in order to e.g. char- 
acterize its atmosphere. For this long-term goal we would like 
to contribute to an "input catalogue" with a significant number 
of entries. Would Doppler spectroscopy have a chance to fulfill 
such an ambitious goal? 

At present, already 2 super-Earths located in the habitable 
zone of their host star have been detected with the HARPS in- 
strument. The first one, G1581 d is part of a multi -planetary sys- 
tem hosting 4 low-mass planets. G1581 d has been detected by 
|Udry et aL| ( 2007 1. As a result of the aliasing with the sideral 
year, its orbital period had to be corrected later (Mayor et al. 
2009a). The minimum mass and orbital period of G1581 d are 
7M e and 66 days. Despite its rather short period the planet ap- 
pears as the first super-Earth discovered in the habitable zone 
due to the very late spectral type and low mass of Gl 581 (M5V, 
0.3 Me). Models of the atmosphere of G1581d have demon- 
strated the possibility of h abitability ([K altenegg er et al.|2011a 



Wor dsworth et al.|[20TT] |Hu & Ding||2011) . The discovery 
of another super-Earth in the habitable zone of the same star 
(G581 g) was claimed by ( Vogtetal.|2010| l. Statistical reanalysis 
of the published velocity data could unfortunately not confirm 
the detection ( |Andrae et al.|2010[ |Gregory|20TTj |Tuomi|201 f 



Doubling the number of available precise HARPS measurements 
|Forveille & al.| ( |2011| ) ruled out the existence of Gl 581 g. 
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A new program to explore the possibility to detect habitable 
planets orbiting solar-type stars with HARPS started 2 years ago 
(Pep e et al.|2Ql 1] >. A measurement strategy is applied to limit the 
influence of stellar noises (acoustic and granulation noises) and 
observations are carried out as often as possible, ideally over sev- 
eral years, in order to extract a potential planetary signal possibly 
hampered by lower frequency noise sources (spots etc). This ob- 
serving strategy is quite demanding in terms of telescope time, 
and forced us to limit the size of this precursor sample to only 
10 solar-type stars. These 10 stars where selected according to 
their proximity (< 16 pes) and for having a low level of chromo- 
spheric activity. After only 2 years, already 5 low-mass planets 
have been discovered orbiting three stars among the ten (Pepe 
et al.||20lT i. We should first mention a 3-planet system hosted 
by the metal deficient ([Fe/H]= -0.40) star G4V HD 20794. For 
these three planets the detected very-low radial-velocity ampli- 
tudes (0.83, 0.56 and 0.85 ms _1 ) are already impressively small 
and correspond to masses of (2.7, 2.4 and 4.8 M ffi ), respectively. 
Another low-mass planet has also been detected ( |Pepe et ah 
20TT) ; it is hosted by the metal-deficient ([Fe/H] = - 0.33) K5V 
star HD85512. Once again, the velocity amplitude is smaller 
than the meter per second (0.77 ± 0.09 ms _1 ). With an orbital pe- 
riod of 58.4 days and a minimum mass of 3.6 M e , HD 855 12 b 
seems to be located inside, but close to the inner boundary, of 
the habitable zone of this K5 star (Kaltene gger et al.|201 lb) . 

Let us notice the proximity to the solar system of G1581 
(6.06 pes) and HD85512 (11.15 pes), and their age very similar 
to the Sun's one (5.8 and 5.6 Gyr). The detection of HD 85512 d 
is close to the HARPS limit of detectability but demonstrates 
the possibility of detecting super-Earths in the habitable zone 
of solar-type stars. To evaluate the sensitivity of Doppler spec- 
troscopy to detect super-Earths in the habitable zone we have 
selected the 10 stars measured with HARPS and having the 
largest number of measurements. When multiple measurements 
are done on the same night, with the purpose to reduce the granu- 
lation noise influence, these measurements are counted for a sin- 
gle data point. All these stars have been measured more than 165 
times during several years. Altogether 29 planets have been dis- 
covered orbiting these ten stars. We have evaluated the detection 
limits for these 1 stars , and sketched them in the m.2 sin i - log P 



plane (see Fig 18 i. Consodering the actual calendar of measure- 



ments of the ongoing HARPS survey, if the number of measure- 
ments is large enough (let say larger than 165 per star), we obtain 
a detection sensitivity close to 100 % for a 10M ffl super-Earth up 
to orbital periods of one year. The corresponding detection prob- 
ability is still close to 20 % for a 3 M ffl planet. 

A new spectrograph called ESPRESSO is presently devel- 
oped for ESO's VLT at the Paranal Observatory. We can also 
mention the development of a northern copy of HARPS to be 
installed on the TNG (3.5-m telescope) at La Palma, Canary 
Islands. This project should contribute to the radial velocity mea- 
surements and mass determination of Kepler mission's planet 
candidates. Part of the HARPS-N observing time will also be 
devoted to the detection of Earth twins orbiting solar-type stars 
in the solar vicinity. We can be confident that in the coming ten 
to twenty years we should have, by the radial-velocity measure- 
ments, an important contribution to a first list of solar-type stars 
with habitable planets and an estimation of their frequency. 

7. Conclusions 

The (still continuing) HARPS survey has already allowed a de- 
tailed study of statistical properties of planetary systems over the 
mass domain from a few Earth masses to several Jupiter masses. 
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Fig. 18. The figure illustrates the limits of detection for very low- 
mass planets based on only the 10 stars of the HARPS sample 
with more than 165 HARPS measurements per star. Mote that 
29 planets have been discovered orbiting these ten stars! 



This survey provides information over three orders of magnitude 
of planetary masses. Thanks to the sub-meter-per-second preci- 
sion of HARPS and the significant number of observing nights 
dedicated to this survey over the last 8 years, exceptional results 
have been obtained. The overwhelming importance of the pop- 
ulation of low mass planets on tight orbit has been quantified. 
About 50% of solar-type stars host planets with mass lower that 
30M ffl on orbits with period shorter than 100 days. The lower 
limit for this frequency value is set by a residual detection bias 
below two to five Earth-masses at periods from few to 100 days. 
Quite obviously any extrapolation to lower masses in that range 
of periods would further increase this occurrence rate. On the 
other hand, because of the possible existence of a population 
of solid planets with rising density below a few Earth-masses 



(Mordasini et al. 2009b) any extrapolation to lower masses is 



highly speculative. 

The time being, the presented mass distribution indicates that 
we are in presence of two populations of planets. They are sep- 
arated by a "gap" in the mass distribution at about 30 Earth- 
masses, which cannot be the consequence of a detection bias. 
These populations have quite different characteristics: 

- The population of gaseous giant planets (GGP) with masses 
above 50 Earth-masses covers the range of periods from days 
to several thousands of days, but the frequency is rising with 
the logarithm of the period log P (Fig.|7]i. These planets ex- 
hibit an upper-mass limit which increases with the orbital 
period (Fig. [6]). The occurrence rate of GGP is about 14% 
(m2sin/ > 50 M e , P < 10 years). Furthermore, the occur- 
rence rate of GGP is strongly correlated with the host star 
metallicity (see black histogram in Fig.[T2]i. In a systematic 
survey to search for planets orbiting metal deficient stars, 
Santos et al.| ( 20lT) have detected only three GGP among 
120 observed stars. All the three GGP have host star metal- 
licities close to -0.5 dex, which was the upper limit set to the 
metallicity of the host stars in that survey. It seems that the 
observations set some lower limit to the host star metallicity 
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(and supposingly to the accretion disk) to allow for the for- 
mation of GGP. Finally, we observe in this sample that the 
distribution of orbital eccentricities exhibits a large scatter 
and may attain values even larger than 0.9. 
- The population of Super-Earth and Neptune-type planets 
(SEN) with masses less than about 30 to 40 Earth-masses be- 
haves quite differently. About 50 percents of solar-type stars 
host at least one super-Earth or Neptune. The mass distri- 
bution is strongly decreasing from 15 to 30M ffl despite the 
fact that they would be detected much more easily (Fig. [8] 



and Fig. 10 1. Also, the SEN population seems to prefer inter- 
mediate orbital periods from 40 to 80 days once corrected 
for (Fig. 12 1, their frequency decreasing on both ends to- 
ward shorter and longer periods. Opposite to GGP, the oc- 
currence rate of SEN does not exhibit a preference for metal 
rich host stars. The difference between both populations of 
planets is striking (see the host-star metallicity histogram 



in Fig. 16 1. The median metallicity for solar-type stars host 
of SEN planets is close to Fe/H = -0.1. A significant num- 
ber of SEN have host stars as metal deficient as -0.4 dex. It 
must be recalled here that only the most massive planet of 
a multi-planetary system is used for that comparative study 
of GGP versus SEN systems. A systematic search is contin- 
uing with HARPS (Pi N.Santos) to explore the occurrence 
of SEN hosted by very metal-deficient stars. At present the 
lower limit of host-star metallicity needed to form a plan- 
etary system remains unknown. A short final word should 
be spent on eccentricities of SEN. They are clearly different 
from zero in many cases, but in contrast to GGP they are 
limited to more modest values up to 0.45. 

It is interesting to note that our estimation for the occurrence 
rate of low-mass planets with periods shorter than 100 days is 
larger than 50 percents. This value is remarkably different from 
the estimation resulting from the Kepler mission. One simple ex- 
planation could come from the fact that Doppler surveys deter- 
mine planetary masses while the transit searches measure plan- 
etary diameters. An interesting alternative to reconcile both es- 



timations of planet occurrence has been proposed by Wolfgang 



& Laughlin (201 1 ): The apparent difference could be the result 
of (at least) two distinct types of planets resulting from multi- 
ple formation mechanisms involved in the production of super- 
Earths. 

The observed planetary populations with their specific char- 
acteristics should help in constraining the physics of planetary 
formation. The discovery, by Doppler spectroscopy, of two plan- 
etary systems with super-Earths in the habitable zone (Gl 58 1 d 
and HD 855 12 b) confirms the potential of that technique to iden- 
tify planets around very close solar-type stars, which are suit- 
able candidates for future experiments designed to detect spec- 
troscopic signatures of life. 
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Table 3. Characteristics of the planetary systems orbiting stars in our sample. 



Name 


Period 


Ecc 


K 


M2sin(i) 


Reference 




[days] 




[ms" 1 ] 


[Earth Mass] 




HD142b 


350.3 


0.25 


33.90 


415.0 


Tinney2002 


HD1461 b 


5.77 


0.00 


2.44 


6.94 


Rivera2010 


HD1461 c 


13.50 


0.00 


1.57 


5.92 


This paper & Segransan201 1 


HD4113b 


526.62 


0.90 


97.09 


523.9 


Tamuz2008 


HD4208 b 


828.0 


0.05 


19.05 


256.6 


Vogt2002 


HD4308 b 


15.62 


0.08 


4.00 


13.45 


Udry2008 


HD6434 b 


21.99 


0.17 


34.20 


126.2 


Mayor2004 


HD7199b 


615 


0.19 


7.8 


92 


Dumusque201 1 


HD7449 b 


1275 


0.82 


41.6 


353 


Dumusque201 1 


HD7449 c 


4046 


0.53 


30 


636 


Dumusque201 1 


HD10180b 


1.17 


0.00 


0.89 


1.52 


Lovis2010 


HD10180c 


5.75 


0.07 


4.53 


13.2 


Lovis2010 


HD10180d 


16.36 


0.16 


2.92 


11.9 


Lovis2010 


HD10180e 


49.75 


0.06 


4.16 


24.8 


Lovis2010 


HD10180f 


122.7 


0.13 


2.94 


23.4 


Lovis2010 


HD10180g 


595 


0.00 


1.62 


22.1 


Lovis2010 


HD10180h 


2149 


0.15 


3.21 


67 


Lovis2010 


HD 10647 b 


1003.0 


0.15 


17.89 


294.0 


Butler2006 


HD11964Ab 


37.94 


0.21 


3.56 


16.8 


Wright2009 


HD11964Ac 


2010 


0.07 


9.74 


176 


Wright2009 


HD 13808 b 


14.18 


0.17 


3.53 


10.33 


This paper & Queloz201 1 


HD 13808 c 


53.83 


0.43 


2.81 


11.55 


This paper & Queloz201 1 


HD16141 b 


75.52 


0.25 


11.98 


79.3 


Marcy2000 


HD16417b 


16.43 


0.22 


3.33 


14.7 


O'Toole2009 


HD 19994 b 


466.2 


0.26 


29.29 


421.7 


Mayor2004 


HD20003 b 


11.84 


0.40 


4.03 


12.00 


This paper & Udry20 11 


HD20003 c 


33.82 


0.16 


2.95 


13.42 


This paper&Udry2011 


HD20781 b 


29.15 


0.11 


3.03 


12.04 


This paper & Udry201 1 


HD20781 c 


85.13 


0.28 


2.88 


15.78 


This paper & Udry201 1 


HD20782 b 


596.2 


0.93 


78.7 


361 


Jones2006 


HD20794 b 


18.31 


0.0 


0.83 


2.7 


Pepe2011 


HD20794 c 


40.11 


0.0 


0.56 


2.4 


Pepe2011 


HD20794 d 


90.30 


0.0 


0.85 


4.8 


Pepe2011 


HD21693 b 


22.65 


0.26 


2.73 


10.22 


This paper & Udry201 1 


HD21693 c 


53.88 


0.24 


4.02 


20.57 


This paper & Udry20 11 


HD23079 b 


730.6 


0.10 


54.90 


776.6 


Tinney2002 


HD27631 b 


2220 


0.17 


27 


540.3 


Marmier2011 


HD28185b 


379.0 


0.05 


163.5 


1842.5 


Santos2001 


HD30562 b 


1157.0 


0.75 


33.70 


423.5 


Fischer2009 


HD3 1527 b 


16.54 


0.13 


3.01 


11.55 


This paper&Udry2011 


HD3 1527 c 


51.28 


0.11 


2.83 


15.82 


This paper & Udry201 1 


HD31527d 


274.49 


0.38 


1.79 


16.50 


This paper & Udry201 1 


HD38858 b 


407.15 


0.27 


2.99 


30.55 


This paper & Benz20 11 


HD39091 b 


2151.0 


0.64 


196.39 


3206.3 


Jones2002 


HD39194b 


5.63 


0.20 


1.95 


3.71 


This paper & Queloz201 1 


HD39194c 


14.02 


0.11 


2.26 


5.94 


This paper & Queloz201 1 


HD39194d 


33.94 


0.20 


1.49 


5.13 


This paper & Queloz201 1 


HD40307 b 


4.31 


0.19 


1.94 


4.07 


Mayor2009 


HD40307 c 


9.62 


0.08 


2.34 


6.51 


Mayor2009 


HD40307 d 


20.43 


0.11 


2.48 


8.85 


Mayor2009 


HD45184b 


5.88 


0.18 


4.09 


11.59 


This paper&Udry2011 


HD45364 b 


226.0 


0.18 


7.44 


61.1 


Correia2009 


HD45364 c 


344.3 


0.12 


22.97 


219 


Correia2009 


[1] Excess residuals. 



Special remarks 



[5] 



[2] Poor coverage or sampling inducing poorly constrained orbital parameters, in particular eccentricity. 

[3] Peak in the RV GLS at periods corresponding to stellar rotation but not overlapping with orbital period of planetary 

companions. Most probably responsible for excess jitter. 
[4] Peak at orbital period in the RV GLS has correspondence with peak at the same period in the R'HK GLS. However, 

after subtraction of the other companions, no correlation is detected between the RV residuals and R'HK. 
[5] Semi-amplitude value of same level than (o-c) rms. 
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Name 


Period 


Ecc 


K 


M2sin(i) 


Reference 


Special remarks 




[days] 




[mr 1 ] 


[Earth Mass] 






HD47186b 


4.08 


0.04 


9.28 


23.0 


Bouchy2009 




HD47186c 


3552 


0.28 


8.12 


185 


Bouchy2009 




HD50499 b 


2457.87 


0.25 


23.02 


554.5 


Vogt2005 




HD5 1608 b 


14.07 


0.15 


4.10 


13.14 


This paper&Udry2011 




HD5 1608 c 


95.41 


0.41 


3.25 


17.97 


This paper & Udry20 11 


[1],[3] 


HD52265 b 


119.29 


0.32 


42.09 


340.5 


Butler2000 




HD65216b 


579 


0.26 


37.8 


449 


Mayor2004 




HD65216 c 


5542 


0.15 


27.5 


712 


Marmier2011 




HD69830 b 


8.66 


0.10 


3.53 


9.86 


Lovis2006 




HD69830 c 


31.58 


0.10 


2.77 


11.98 


Lovis2006 




HD69830 d 


199.56 


0.16 


2.12 


16.73 


Lovis2006 




HD70642 b 


2068.0 


0.03 


30.39 


606.9 


Carter2003 




HD75289 b 


3.50 


0.03 


54.90 


146.2 


Udry2000 




HD82943 c 


219.5 


0.35 


66.0 


632.2 


Mayor2004 




HD82943 b 


441.2 


0.21 


43.59 


551.0 


Mayor2004 




HD83443 b 


2.98 


0.01 


56.20 


125.8 


Butler2002 




HD85390 b 


809.4 


0.44 


3.95 


43.1 


Mordasini2010 




HD85512b 


58.43 


0.11 


0.76 


3.6 


Pepe2011 




HD86226 b 


1681 


0.11 


17.2 


322 


Arriagada2010 




HD90156 b 


49.88 


0.46 


3.66 


16.7 


Mordasini2010 




HD92788 b 


325.81 


0.33 


106.0 


1132.7 


Fischer2001 




HD93083 b 


144.2 


0.07 


19.9 


136.5 


Lovis2005 




HD93385 b 


13.18 


0.15 


2.21 


8.36 


This paper & Queloz201 1 




HD93385 c 


46.02 


0.24 


1.82 


10.12 


This paper & Queloz201 1 




HD96700 b 


8.12 


0.10 


3.02 


9.02 


This paper & Queloz201 1 


[1],[3] 


HD96700 c 


103.49 


0.37 


1.98 


12.76 


This paper & Queloz201 1 


[1],[3] 


HD98649 b 


10400 


0.86 


135 


2225 


Marmier2011 




HD101930b 


70.49 


0.07 


19.1 


104.2 


Lovis2005 




HD102117b 


20.82 


0.07 


11.1 


48.0 


Lovis2005 




HD 102365 b 


122.1 


0.34 


2.40 


16.2 


Tinney201 1 




HD 104067 b 


55.83 


0.06 


10.9 


51.4 


Segransan201 1 




HD106515Ab 


3630 


0.60 


174 


3337 


Marmier2011 




HD108147b 


10.89 


0.52 


25.10 


82.0 


Pepe2002 




HD1 11232 b 


1118.0 


0.18 


162.0 


2199.7 


Mayor2004 




HD1 14386 b 


445 


0.12 


11.9 


119 


Mayor2004 




HD 114386 c 


1046 


0.06 


28.1 


377 


Marmier2011 




HD1 14729 b 


1114.0 


0.16 


18.79 


300.3 


Butler2003 




HD1 14783 b 


493.7 


0.14 


31.89 


351.2 


Vogt2002 




HD115617b 


4.21 


0.20 


2.59 


5.83 


Vogt2010 




HD1 15617 c 


38.07 


0.18 


4.68 


22.20 


Vogt2010 




HD115617d 


123.01 


0.34 


3.25 


22.8 


Vogt2010 




HD 117207 b 


2597.0 


0.14 


26.60 


578.1 


Marcy2005 




HD1 17618 b 


25.82 


0.41 


12.80 


56.1 


Tinney2005 




HD121504b 


63.33 


0.02 


55.79 


388.5 


Mayor2004 




HD 126525 b 


948.12 


0.13 


5.11 


71.33 


This paper & Benz201 1 




HD 134060 b 


3.27 


0.40 


4.97 


11.17 


This paper & Udry20 11 




HD 134060 c 


1160.97 


0.75 


4.06 


47.87 


This paper&Udry2011 




HD 134606 b 


12.08 


0.15 


2.68 


9.27 


This paper & Queloz201 1 




HD 134606 c 


59.51 


0.29 


2.17 


12.14 


This paper & Queloz201 1 


[4] 


HD 1 34606 d 


459.26 


0.46 


3.66 


38.52 


This paper & Queloz201 1 




HD 134987 b 


258.18 


0.23 


49.5 


496.9 


Vogt2000 




HD 134987 c 


5000.0 


0.11 


9.30 


255.8 


Vogt2000 




HD 136352 b 


11.57 


0.18 


1.77 


5.28 


This paper & Udry20 11 


[5] 


[1] Excess residuals. 



[2] Poor coverage or sampling inducing poorly constrained orbital parameters, in particular eccentricity. 

[3] Peak in the RV GLS at periods corresponding to stellar rotation but not overlapping with orbital period of planetary 

companions. Most probably responsible for excess jitter. 
[4] Peak at orbital period in the RV GLS has correspondence with peak at the same period in the R'HK GLS. However, 

after subtraction of the other companions, no correlation is detected between the RV residuals and R'HK. 
[5] Semi-amplitude value of same level than (o-c) rms. 
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Table 3. continued. 



Name 


Period 


Ecc 


K 


M2sin(i) 


Reference Special remarks 




[days] 




[ms ] 


[Earth Mass] 




HD 136352 c 


27.58 


0.16 


2.82 


11.38 


This paper & Udry201 1 [5] 


HD 136352 d 


106.72 


0.43 


1.68 


9.58 


This paper & Udry201 1 [5] 


HD137388 b 


330 


0.36 


7.9 


71 


Dumusque201 1 


HD141937b 


653.21 


0.40 


234.5 


3011.5 


Udry2002 


HD 142022 A b 


1928.0 


0.52 


92.0 


1419.9 


Eggenberger2006 


HD147018 b 


44.23 


0.46 


145.33 


676.1 


Segransan2010 


HD147018 c 


1008.0 


0.13 


141.19 


2095.6 


Segransan2010 


HD150433 b 


1096.27 


-0.00 


3.85 


53.48 


This paper & Benz201 1 [2] 


HD 15408 8 b 


18.59 


0.38 


1.78 


6.14 


This paper & Queloz201 1 


HD 156846 b 


359.51 


0.84 


464.29 


3498.5 


Tamuz2008 


HD157172b 


105.0 


0.33 


5.46 


34.6 


This paper & Benz201 1 


HD 160691 b 


9.64 


0.12 


3.05 


10.6 


Pepe2007 


HD 160691 c 


313.2 


0.04 


16.9 


189.1 


Pepe2007 


HD 160691 d 


648.7 


0.18 


39.0 


546.9 


Pepe2007 


HD 160691 e 


8723 


0.43 


25.8 


790 


Pepe2007 


HD 166724 b 


8100 


0.77 


72 


1310 


Marmier2011 


HD 168443 b 


58.11 


0.52 


475.53 


2465.3 


Marcyl999 


HD 168443 c 


1748.15 


0.21 


298.13 


5571.5 


Marcy2001 


HD 168746 b 


6.40 


0.10 


28.60 


77.9 


Pepe2002 


HD 169830 b 


225.62 


0.31 


80.69 


918.4 


Naef2001 


HD169830 c 


2100.0 


0.33 


54.29 


1291.3 


Mayor2004 


HD 179949 b 


3.09 


0.02 


112.59 


286.7 


Tinney2001 


HD181433 b 


9.37 


0.42 


2.89 


7.32 


Bouchy2009 


HD 18 1433 c 


1019 


0.25 


17.2 


222 


Bouchy2009 


HD 18 1433 d 


3201 


0.11 


9.49 


184 


Bouchy2009 


HD 187085 b 


986.0 


0.46 


17.0 


255.4 


Jones2006 


HD 189567 b 


14.27 


0.23 


3.02 


10.03 


This paper & Queloz201 1 


HD192310b 


74.72 


0.13 


3.00 


16.9 


Howard2010 


HD192310c 


525.8 


0.32 


2.27 


24 


Pepe2011 


HD 196050 b 


1378.0 


0.22 


49.70 


903.7 


Jones2002 


HD 196067 b 


4100 


0.63 


112 


2257 


Marmier2011 


HD204313 c 


34.88 


0.13 


3.28 


16.9 


This paper & Segransan201 1 [4] 


HD204313 b 


2132 


0.11 


72.9 


1479 


Segransan2009 


HD204941 b 


1733 


0.37 


5.94 


85 


Dumusque201 1 


HD208487 b 


130.08 


0.23 


19.70 


162.8 


Tinney2005 


HD2 10277 b 


442.19 


0.47 


38.93 


404.5 


Marcyl999 


HD2 13240 b 


882.7 


0.42 


96.59 


1440.6 


Santos2001 


HD215152b 


7.28 


0.34 


1.26 


2.77 


This paper & Queloz201 1 [5] 


HD215152c 


10.86 


0.38 


1.26 


3.09 


This paper & Queloz201 1 [5] 


HD2 15456 b 


192.0 


0.13 


3.62 


32.6 


This paper & Benz201 1 


HD2 15456 c 


2268 


0.17 


3.77 


76.8 


This paper & Benz201 1 


HD2 16435 b 


1311.0 


0.07 


19.60 


386.1 


Jones2003 


HD2 16437 b 


1 1 CI A 

1353.0 


0.31 


OA A 

39.0 


CO A 1 

689.1 


1* „ r\ A AO 

Jones2002 


HD2 16770 b 


118.45 


0.37 


30.89 


205.6 


Mayor2004 


HD217107b 


7.12 


0.12 


139.20 


445.3 


Fischerl999 


HD217107c 


4270.0 


0.51 


35.70 


831.3 


Vogt2005 


HD2 18566 b 


225.7 


0.30 


8.30 


67.7 


Meschiari201 1 


HD220689 b 


2191 


0.2 


19. 


377 


Marmier2011 


HD222582 b 


572.38 


0.72 


276.29 


2425.1 


Vogt2000 



[1] Excess residuals. 

[2] Poor coverage or sampling inducing poorly constrained orbital parameters, in particular eccentricity. 

[3] Peak in the RV GLS at periods corresponding to stellar rotation but not overlapping with orbital period of planetary 



companions. Most probably responsible for excess jitter. 
[4] Peak at orbital period in the RV GLS has correspondence with peak at the same period in the R'HK GLS. However, 

after subtraction of the other companions, no correlation is detected between the RV residuals and R'HK. 
[5] Semi-amplitude value of same level than (o-c) rms. 
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M. Mayor et al.: The HARPS search for southern extra-solar planets 



Table 4. Orbital and physical parameters of the planets orbiting HD 1461 as obtained from a Keplerian fit to the data. 



Parameter 


[unit] 


HD 1461 c 


V 


[kms 4 ] 


-10.0670 (±0.0002) 


D 

r 


[days] 


1 J. JUJ 






/_i_n nn/i\ 


€ 




n n 






(fixed) 


K 


[ms- 1 ] 


1.57 






(±0.19) 


111 CI Tl 1 
III hill i 


\M 1 

t™ffij 








(±0.76) 


a 


[AU] 


0.1117 






(±0.0018) 


N 

1 * meas 




167 


Span 


[days] 


2856 


rms 


[ms- 1 ] 


1.70 


xl 




3.61 



Table 5. Orbital and physical parameters of the planets orbiting HD 13808 as obtained from a Keplerian fit to the data. 



Parameter 


[unit] 


HD 13808 b HD 13808 c 


V 


[kms- 1 ] 


41.0949 (±0.0005) 


P 


[days] 


14.182 


53.83 






(±0.005) 


(±0.11) 


e 




0.17 


0.43 






(±0.07) 


(±0.20) 


K 


[ms-'] 


3.53 


2.81 






(±0.29) 


(±0.46) 


m sin i 


[M e ] 


10.33 


11.55 






(±0.92) 


(±1.62) 


a 


[AU] 


0.1017 


0.2476 






(±0.0016) 


(±0.0041) 


^meas 




133 




Span 


[days] 


2770 




rms 


[ms" 1 ] 


1.93 





410 



Table 6. Orbital and physical parameters of the planets orbiting HD 20003 as obtained from a Keplerian fit to the data. 



Parameter 


[unit] 


HD 20003 b HD 20003 c 


V 


[kms" 1 ] 


-16.1040 (± 


0.0006) 


P 


[days] 


11.849 


33.823 






(±0.003) 


(±0.065) 


e 




0.40 


0.16 






(±0.08) 


(±0.09) 


K 


[ms- 1 ] 


4.03 


2.95 






(±0.33) 


(±0.28) 


m sin i 


[M e ] 


12.00 


13.42 






(±0.97) 


(±1.28) 


a 


[AU] 


0.0974 


0.1961 






(±0.0016) 


(±0.0032) 


N 

iv meas 




104 




Span 


[days] 


2770 




rms 


[ms- 1 ] 


1.56 
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M. Mayor et al.: The HARPS search for southern extra-solar planets 



Table 7. Orbital and physical parameters of the planets orbiting HD 2078 1 as obtained from a Keplerian fit to the data. 



Parameter 


[unit] 


HD 20781 b HD 20781 c 


V 


[kms- 1 ] 


40.3671 (±0.0005) 


D 

r 


[days] 


29.15 


85.13 






(±0.02) 


(±0.12) 


€ 




0.11 


0.28 






(±0.06) 


(±0.09) 


K 


[ms- 1 ] 


3.03 


2.88 






(±0.26) 


(±0.23) 


m sin i 


[Mb] 


12.04 


15.78 






(±1.12) 


(±1.21) 


a 


[AU] 


0.1690 


0.3456 






(±0.0028) 


(±0.0057) 


Apneas 




96 




Span 


[days] 


2647 




rms 


[ms" 1 ] 


1.12 





2d L39 



Table 8. Orbital and physical parameters of the planets orbiting HD 21693 as obtained from a Keplerian fit to the data. 



Parameter 


[unit] 


HD 21693 b HD 21693 c 


V 


[kms- 1 ] 


39.7685 (±0.0008) 


P 


[days] 


22.656 


53.88 






(±0.024) 


(±0.07) 


e 




0.26 


0.24 






(±0.17) 


(±0.09) 


K 


[ms-'] 


2.73 


4.02 






(±0.33) 


(±0.35) 


m sin i 


[M e ] 


10.22 


20.57 






(±1.46) 


(±1.80) 


a 


[AU] 


0.1484 


0.2644 






(±0.0024) 


(±0.0044) 


Apneas 




128 




Span 


[days] 


2813 




rms 


[ms" 1 ] 


2.02 





xl 5,16 



Table 9. Orbital and physical parameters of the planets orbiting HD 31527 as obtained from a Keplerian fit to the data. 



Parameter 


[unit] 


HD 31527 b 


HD 31527 c 


HD 31527 d 


V 


[kms" 1 ] 




25.7391 (±0.0004) 




P 


[days] 


16.546 


51.28 


274.5 






(±0.007) 


(±0.09) 


(±7.8) 


e 




0.13 


0.11 


0.38 






(±0.05) 


(±0.07) 


(±0.25) 


K 


[ms- 1 ] 


3.01 


2.83 


1.79 






(±0.18) 


(±0.17) 


(±0.68) 


m sin ( 


[M ffi ] 


11.55 


15.82 


16.50 






(±0.79) 


(±1.10) 


(±3.04) 


a 


[AU] 


0.1253 


0.2665 


0.818 






(±0.0020) 


(±0.0044) 


(±0.020) 


N 

1 v meas 






167 




Span 


[days] 




2719 




rms 


[ms- 1 ] 




1.35 




x] 






2.14 





M. Mayor et al.: The HARPS search for southern extra-solar planets 
Table 10. Orbital and physical parameters of the planets orbiting HD 38858 as obtained from a Keplerian fit to the data. 



Parameter 


[unit] 


HD 38858 b 


V 


[kms- 1 ] 


31.6383 (±0.0004) 


D 

r 


[days] 


A(\l 1 






(±4.3; 


€ 




n ?7 






(±0.17) 


K 


[ms- 1 ] 


2.99 






(±0.33) 


til alll i 


\M 1 

t™eJ 








(±4.11) 


a 


[AU] 


1.038 






(±0.019) 


N 

1 ' meas 




52 


Span 


[days] 


2964 


rms 


[ms- 1 ] 


1.38 


xl 




2.53 



Table 11. Orbital and physical parameters of the planets orbiting HD 39194 as obtained from a Keplerian fit to the data. 



Parameter 


[unit] 


HD 39194 b 


HD 39194 c 


HD 39194 d 


V 


[kms- 1 ] 




14.1692 (±0.0005) 




P 


[days] 


5.6363 


14.025 


33.941 






(±0.0008) 


(±0.005) 


(±0.035) 


e 




0.20 


0.11 


0.20 






(±0.10) 


(±0.06) 


(±0.16) 


K 


[ms- 1 ] 


1.95 


2.26 


1.49 






(±0.16) 


(±0.15) 


(±0.17) 


m sin ( 


[MJ 


3.72 


5.94 


5.14 






(±0.33) 


(±0.47) 


(±0.66) 


a 


[AU] 


0.0519 


0.0954 


0.1720 






(±0.0008) 


(±0.0016) 


(±0.0029) 








133 




Span 


[days] 




2717 




rms 


[ms- 1 ] 




1.11 




xl 






1.30 





Table 12. Orbital and physical parameters of the planets orbiting HD45184 as obtained from a Keplerian fit to the data. 



Parameter 


[unit] 


HD 45184 b 


V 


[kms" 1 ] 


-3.7584 (±0.0008) 


P 


[days] 


5.8872 






(±0.0015) 


e 




0.30 






(±0.19) 


K 


[ms- 1 ] 


4.77 






(±1.18) 


m sin ( 


[M ffi ] 


12.73 






(±1.67) 


a 


[AU] 


0.0638 






(±0.0010) 


N 

1 v meas 




82 


Span 


[days] 


2738 


rms 


[ms- 1 ] 


3.81 


xl 




21.86 
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M. Mayor et al.: The HARPS search for southern extra-solar planets 



Table 13. Orbital and physical parameters of the planets orbiting HD 51608 as obtained from a Keplerian fit to the data. 



Parameter 


[unit] 


HD 51608 b HD 51608 c 


V 


[kms -1 ] 


39.9773 (±0.0008) 


D 

r 


[days] 


14.070 


95.42 






(±0.004) 


(±0.39) 


€ 




0.15 


0.41 






(±0.06) 


(±0.18) 


K 


[ms- 1 ] 


4.10 


3.25 






(±0.27) 


(±0.61) 


m sin i 


[M e ] 


13.14 


17.97 






(±0.98) 


(±2.61) 


a 


[AU] 


0.1059 


0.379 






(±0.0017) 


(±0.006) 


Apneas 




118 




Span 


[days] 


2697 




rms 


[ms" 1 ] 


1.82 





2d 3.95 



Table 14. Orbital and physical parameters of the planets orbiting HD 93385 as obtained from a Keplerian fit to the data. 



Parameter 


[unit] 


HD 93385 b HD 93385 c 


V 


[kms- 1 ] 


47.5791 (±0.0005) 


P 


[days] 


13.186 


46.025 






(±0.006) 


(±0.073) 


e 




0.15 


0.24 






(±0.11) 


(±0.18) 


K 


[ms-'] 


2.21 


1.82 






(±0.23) 


(±0.96) 


m sin i 


[M e ] 


8.36 


10.12 






(±0.88) 


(±1.47) 


a 


[AU] 


0.1116 


0.2570 






(±0.0018) 


(±0.0043) 


Apneas 




127 




Span 


[days] 


2733 




rms 


[ms" 1 ] 


1.42 





xl 238 



Table 15. Orbital and physical parameters of the planets orbiting HD 96700 as obtained from a Keplerian fit to the data. 



Parameter 


[unit] 


HD 96700 b HD 96700 c 


V 


[kms" 1 ] 


12.8584 (±0.0071) 


P 


[days] 


8.1256 


103.49 






(±0.0013) 


(±0.58) 


e 




0.10 


0.37 






(±0.05) 


(±0.19) 


K 


[ms- 1 ] 


3.02 


1.98 






(±0.18) 


(±0.37) 


m sin i 


[M ffi ] 


9.03 


12.76 






(±0.63) 


(±1.63) 


a 


[AU] 


0.0774 


0.422 






(±0.0012) 


(±0.007) 


N 

iv meas 




146 




Span 


[days] 


2717 




rms 


[ms- 1 ] 


1.48 




xl 




2.93 
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Table 16. Orbital and physical parameters of the planets orbiting HD 126525 as obtained from a Keplerian fit to the data. 



Parameter 


[unit] 


HD 126525 b 


V 


[kms^ 1 ] 


13.1384 (±0.0006) 


D 

r 


[days] 


Q/1 S 1 






(+ZZ.U) 


€ 




1 ^ 






(±0.07) 


K 


[ms- 1 ] 


5.11 






(±0.34) 


til alll i 


t™el 


/ 1 .j j 






(±5.78) 


a 


[AU] 


1.811 






(±0.041) 


N 

1 ' meas 




45 


Span 


[days] 


2575 


rms 


[ms- 1 ] 


1.27 


xl 




2.10 



Table 17. Orbital and physical parameters of the planets orbiting HD 134060 as obtained from a Keplerian fit to the data. 



Parameter 


[unit] 


HD 134060 b HD 134060 c 


V 


[kms- 1 ] 


37.9893 (±0.0026) 


P 


[days] 


3.2700 




1161 






(±0.0002) 




(±27) 


e 




0.40 




0.75 






(±0.04) 




(±0.19) 


K 


[ms- 1 ] 


4.97 




4.1 






(±0.23) 




(±1.8) 


m sin ( 


[MJ 


11.17 




47.9 






(±0.66) 




(±22.5) 


a 


[AU] 


0.0444 




2.226 






(±0.0007) 




(±0.051) 


N meas 






100 




Span 


[days] 




2615 




rms 


[ms" 1 ] 




1.39 










2.65 





Table 18. Orbital and physical parameters of the planets orbiting HD 134606 as obtained from a Keplerian fit to the data. 



Parameter 


[unit] 


HD 134606 b 


HD 134606 c 


HD 134606 d 


V 


[kms" 1 ] 




2.0510 (±0.0006) 




P 


[days] 


12.083 


59.52 


459.3 






(±0.010) 


(±0.17) 


(±8.3) 


e 




0.15 


0.29 


0.46 






(±0.10) 


(±0.20) 


(±0.09) 


K 


[ms- 1 ] 


2.68 


2.17 


3.66 






(±0.25) 


(±0.35) 


(±0.62) 


m sin i 


[M e ] 


9.27 


12.14 


38.52 






(±0.95) 


(±1.67) 


(±4.12) 


a 


[AU] 


0.1023 


0.2962 


1.157 






(±0.0017) 


(±0.0049) 


(±0.024) 


N 

1 v meas 






113 




Span 


[days] 




2548 




rms 


[ms- 1 ] 




1.66 




xl 






3.96 
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Table 19. Orbital and physical parameters of the planets orbiting HD 136352 as obtained from a Keplerian fit to the data. 



Parameter 


[unit] 


HD 136352 b 


HD 136352 c 


HD 136352 d 


V 


[kms- 1 ] 




-68.7120 (±0.0005) 




D 

r 


[days] 


1 1 <T7 


in ^co 

z / joz 


1UD. 1 L 






(+U.UUo) 


/ i A AT3\ 

(+U.U2J) 


(±1.03) 


B 




U.lO 


U.1U 








(±0.14) 


(±0.07) 


(±0.24) 


K 


[ms- 1 ] 


1.77 


2.82 


1.68 






(±0.22) 


(±0.23) 


(±0.47) 


ill c in 7 
ill sill l 






1 1 .JO 


y — >y 






(±0.62) 


(±0.10) 


(±1.86) 


a 


[AU] 


0.0933 


0.1665 


0.411 






(±0.0015) 


(±0.0028) 


(±0.007) 


N 

1 ' meas 






129 




Span 


[days] 




2601 




rms 


[ms- 1 ] 




1.35 




xl 






2.72 





Table 20. Orbital and physical parameters of the planets orbiting HD 150433 as obtained from a Keplerian fit to the data. 



Parameter 


[unit] 


HD 150433 b 


V 


[kms- 1 ] 


-40.1141 (±0.0003) 


P 


[days] 


1096 






(±27) 


e 




0.0 






(fixed) 


K 


[ms- 1 ] 


3.85 






(±0.42) 


m sin i 


[M m ] 


53.5 






(±6.3) 


a 


[AU] 


1.930 






(±0.045) 


N meas 




51 


Span 


[days] 


2129 


rms 


[ms- 1 ] 


1.73 


xl 




3.75 



Table 21. Orbital and physical parameters of the planets orbiting HD 154088 as obtained from a Keplerian fit to the data. 



Parameter 


[unit] 


HD 154088 b 


V 


[kms- 1 ] 


14.2972 (±0.0003) 


P 


[days] 


18.596 






(±0.021) 


e 




0.38 






(±0.15) 


K 


[ms- 1 ] 


1.78 






(±0.31) 


m sin ( 


[M e ] 


6.15 






(±0.86) 


a 


[AU] 


0.1316 






(±0.0021) 


N 

1 ' meas 




112 


Span 


[days] 


1924 


rms 


[ms- 1 ] 


1.24 


xl 




2.21 



M. Mayor et al.: The HARPS search for southern extra-solar planets 
Table 22. Orbital and physical parameters of the planets orbiting HD 157172 as obtained from a Keplerian fit to the data. 



Parameter 


[unit] 


HD 157172 b 


V 


[kms 4 ] 


-78.9210 (±0.0006) 


D 

r 


[days] 


1U4.84 






(±0.13) 


€ 










(±0.05) 


K 


[ms- 1 ] 


6.42 






(±0.43) 


111 CI Tl 1 
III hill i 


\M 1 

t™ffij 


38 1 

JO. 1 






(±2.6) 


a 


[AU] 


0.416 






(±0.007) 


N 

1 * meas 




82 


Span 


[days] 


2157 


rms 


[ms- 1 ] 


1.90 


xl 




13.03 



Table 23. Orbital and physical parameters of the planets orbiting HD 189567 as obtained from a Keplerian fit to the data. 



Parameter 


[unit] 


HD 189567 b 


V 


[kms- 1 ] 


-10.4792 (±0.0004) 


P 


[days] 


14.275 






(±0.005) 


e 




0.23 






(±0.14) 


K 


[ms- 1 ] 


3.02 






(±0.33) 


m sin i 


[M m ] 


10.03 






(±1.04) 


a 


[AU] 


0.1099 






(±0.0018) 


N meas 




166 


Span 


[days] 


2818 


rms 


[ms- 1 ] 


2.64 


xl 




9.09 



Table 24. Orbital and physical parameters of the planets orbiting HD 204313 as obtained from a Keplerian fit to the data. 



Parameter 


[unit] 


HD 204313 c 


V 


[kms" 1 ] 


-9.7419 (±0.0021) 


P 


[days] 


34.873 






(±0.039) 


e 




0.17 






(±0.09) 


K 


[ms- 1 ] 


3.36 






(±0.35) 


m sin ( 


[M ffi ] 


17.15 






(±1.71) 


a 


[AU] 


0.2103 






(±0.0035) 


N 

1 v meas 




67 


Span 


[days] 


1547 


rms 


[ms- 1 ] 


1.08 


xl 




1.78 
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Table 25. Orbital and physical parameters of the planets orbiting HD 215152 as obtained from a Keplerian fit to the data. 



Parameter 


[unit] 


HD 215152 b HD 215152 c 


V 


[kms- 1 ] 


-13.6810(±0.0012) 


P 

r 


[days] 


7.283 




10.866 






(±0.006) 




(±0.014) 


€ 




0.34 




0.38 






(±0.27) 




(±0.23) 


K 


[ms- 1 ] 


1.26 




1.26 






(±0.36) 




(±0.32) 


11 1 c 1 Tl 1 

til alii i 


\M 1 

t™ffij 


2.78 




3.10 






(±0.47) 




(±0.48) 


a 


[AU] 


0.0652 




0.0852 






(±0.0010) 




(±0.0014) 


N 






171 




Span 


[days] 




2927 




rms 


[ms- 1 ] 




1.33 




xl 






2.15 





Table 26. Orbital and physical parameters of the planets orbiting HD 215456 as obtained from a Keplerian fit to the data. 



Parameter 


[unit] 


HD 215456 b HD 215456 c 


V 


[kms- 1 ] 


-18.8541 (±0.0006) 


P 


[days] 


191.99 


2277 






(±0.73) 


(±67) 


e 




0.15 


0.19 






(±0.10) 


(±0.11) 


K 


[ms- 1 ] 


3.62 


3.89 






(±0.31) 


(±0.44) 


m sin ( 


[MJ 


32.21 


78.37 






(±2.92) 


(±8.86) 


a 


[AU] 


0.652 


3.394 






(±0.011) 


(±0.088) 


N meas 




96 




Span 


[days] 


2615 




rms 


[ms- 1 ] 


1.79 





